(N 



< 



Thermal abundance of non-relativistic relics with 
Sommerfeld enhancement 



a -43, Mitsuru Kakizaki c 



Hoernisa Iminniyaz a ' Mitsuru Kakizaki 

a Center for High Energy Physics, Peking University, Peking, 100871, China 
6 School of Physics Science and Technology, Xinjiang University, Urumqi, 830046, 

China 



C LAPTH, Universite de Savoie, CNRS, B.P. 110, F-74941 Annecy-le-Vieux Cedex, 



France 



o. 

, Abstract 

ji ■ 

O !■ We propose an analytic treatment for computing the relic abundances of 

O ■ non-relativistic particles whose annihilation rate at chemical decoupling is in- 

creased by Sommerfeld enhancement. We find approximate rational functions 
that closely fit the thermal average of Sommerfeld-enhanced cross sections in 
the massless limit of force carriers for s- and p-wave annihilations. We demon- 
strate that, with the approximate thermally-averaged cross sections imple- 
mented, the standard analytic method for the final relic abundances provides 
£S| ' accuracy to within 1% even for the case of Sommerfeld enhancement. 
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1 Introduction 



The determination of the relic abundance of particles which decouple from thermal 
equilibrium in the early universe is indispensable for understanding the history of 
the universe. Important examples include the computation of the cosmological dark 
matter abundance, which provides us a crucial hint for screening dark matter can- 
didates and models beyond the standard model (SM) of particle physics, as well 
as cosmological scenarios in the early universe. Among many particles proposed, 
stable or long-lived weakly interacting massive particles (WIMPs) with weak-scale 
mass are excellent candidates because, assuming the thermal production scenario, 
the predicted relic density coincides with the dark matter density [H [2] . The value 
extracted from the Wilkinson Microwave Anisotropy Probe (WMAP) data is [3], 



where Odm is the ratio of the dark matter mass density to the critical density, and 
h = 0.710 d= 0.025 is the scaled Hubble parameter in units of 100 km sec" 1 Mpc" 1 . 
The uncertainty will be improved by the PLANCK satellite [3]. 

In order to calculate the number density of relic particles x accurately, in prin- 
ciple one must solve the Boltzmann equation, which describes the evolution of the 
distribution function. In the standard thermal WIMP production scenario, where 
decoupling occurs in the radiation-dominated epoch, the particle number density 
is determined only by the thermal average of the product of the annihilation cross 
section a and the relative velocity of the annihilating particles v. In many cases, 
the thermally-averaged annihilation cross section can be expanded in a power se- 
ries: (cry) = a + 6b/ x + • • • , where x is the ratio of the particle's mass m x to the 
temperature T, leading to simple analytic formulas for the final abundance [5j |6], 
although there are some exceptional cases [7] . The desired cross section for reconcil- 
ing with the WMAP range is found to be approximately (av) ~ 3 x 10~ 26 cm 3 sec -1 
for the temperature at which WIMPs decouple from the thermal bath. Analytic 
methods for the abundance of relic particles have been developed also in various 
non-standard cosmological scenarios where the relic abundance is increased or de- 
creased due to a low reheat temperature, the late decay of a scalar field, entropy 
production at late times, modification of the Hubble expansion rate, or their com- 
bination [H [91 [101 HU [Til [131 HH ISl ESI IZl [HI lU ED]. Such analytic approaches 
enable us to estimate the relic abundance without tedious numerical computations. 




0.1109 ±0.0056 




1 



Recently, it was pointed out that the annihilation rate of dark matter particles can 
be significantly altered by the so-called Sommerfeld enhancement at low velocities 
[21| |22| [23j 121]. Anomalous excesses of cosmic positrons reported by PAMELA 
[25] , ATIC [26J and FERMI [21] have motivated us to investigate the Sommerfeld 
effect because the resulting annihilation cross section is extremely boosted for lower 
velocities [28 ] 129 ] 130] . However, it should be also emphasized that the relic abundance 
can be significantly reduced by the Sommerfeld effect at chemical decoupling [23] 
[21]. There is also a possibility of a change in the relic abundance after kinetic 
decoupling due to chemical recoupling of the annihilation interactions |31j. (Semi- 
analytic treatment for relic abundances for the case of Sommerfeld-enhanced s-wave 
annihilations has been discussed [321 E21 EH US] . It was shown that to a very good 
approximation the standard analytic method reproduces numerically computed relic 
abundances for Sommerfeld-enhanced s-wave annihilations in the limit a/v 3> 1, 
where a is a coupling constant between the WIMP and force carrier [33]. However 
until now, no analytic formula for approximating relic abundances has been discussed 
for arbitrary a/v, even for the case of massless force carriers. 

In this paper, we address the relic abundance of non-relativistic particles whose 
annihilation rate is altered by the Sommerfeld enhancement when the relic parti- 
cles decouple from the thermal background!. We find highly accurate functions that 
describes the transition from non-enhanced thermally averaged cross sections to 1/v- 
enhancement not only for s-wave but also for p-wave annihilations. We then show 
that the standard analytic method for the final relic abundances provides accuracy 
to within 1% even in the range where 1/v approximation does not work. The de- 
rived approximate formula is a powerful tool for estimating relic abundances before 
potential chemical recoupling. 

This paper is arranged as follows. In Section 2, we discuss our method for approx- 
imating the thermal average of Sommerfeld-enhanced annihilation cross sections. In 
Section 3, we describe the standard method for deriving the relic density including 
Sommerfeld enhancement, and compare our analytic results to numerical computa- 
tions. Section 4 is devoted to our conclusions. 



f In Ref. EH the case with massive force carriers is discussed. 
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2 Thermally-averaged annihilation cross section 



In this section, we discuss approximated expressions for the thermal average of 
Sommerfeld-enhanced WIMP annihilation cross sections. 

When WIMPs decouple from thermal background, they are non-relativistic. In 
the absence of force carrier, the annihilation cross section of WIMPs can be expanded 
with respect to the relative velocity v, 

a v = a + bv 2 + 0(v 4 ) , (2) 

where a, b are constants. For s-wave annihilation, a gives the dominant contribution 
to the annihilation of WIMPs. If the s-wave contribution is suppressed, b is described 
by the p-wave contribution. When massless force carriers mediate interactions be- 
tween annihilating particles, the annihilation cross section is enhanced by the factor 
[36], 

„2" 



s, 



n * 2 +, 

s=l x 



a 

~~2 



e ™/v na / V (Q) 

smh{7Ta/v){l\) 2 ' y ' 



for /-partial wave. Here, a is a coupling strength. For example, for the case where the 
annihilation rate of a fermion-antifermion pair is enhanced by scalar boson exchanges, 
the coupling strength is given by a = f 2 /(47t), with / being the Yukawa coupling 
constant. In this paper, we focus on the first two lowest modes, 

2ira/v _ f a?\ 2na/v 

s 1 _ e -2na/v ' °P \ ^ v 2 J l_ e -2ira/v " W 

The case of a/v <C 1 results in the usual non-enhanced annihilation cross sections as 
Si approaches to unity. In the opposite limit a/v ^> 1, the Sommerfeld enhancement 
factor is simplified down to S s = luajv for s-wave annihilation, and to S p = 2ira 3 /v 3 
for p-wave annihilation. At the leading order of the v 2 expansion, we can parametrize 
the Sommerfeld-enhanced annihilation cross sections as 

\ aS s (s-wave) 
<jv = < (5) 
I bv 2 S p (p-wave) . 

Inclusion of higher order terms is a trivial extension of the procedure we will present. 

Let us consider the enhancement of the WIMP annihilation by introducing the 
boost factor B = (av) / (ctqv) . The thermal average of the s-wave annihilation cross 
section is given by 

(av) = a(S s) = a ^ Hdv v* - ^ , . (6) 

\ / n s ' 2 Jtt L 1 - e- 2 ™/^ v ; 



Introducing further the variables y = a^pnx and t = v v / x/(2 A /7r), the boost factor 
is solely described by y, 



POD 

B s (y) = (S a ) = Any / dt t 
Jo 



Similarly, for p-wave annihilation, we obtain 



Q 3 J V 47r ^ 2 / 1 - e ~ v/t ' 

In the case where the coupling constant a is small or the temperature is high 
enough to suppress Sommerfeld enhancement, we can expand the boost factor in the 
Taylor series, 

£ T aylor,s(y) = 1 + V + ^V* > 

2 3 + 7T 2 „ 

Bt^M = l + -y + ——y 2 . (9) 

O l07T 

In the opposite limit, where the cross section is enhanced by 1/v, e~ y ^ in the de- 
nominators of Eqs.([7D and (JSJ) are negligible, leading to 

B 1/VjS (y) = 2y , B 1/v>p (y) = t- + ±y. (10) 

As a simple function that connects the two limiting results for y — y and y — > oo, 
we propose the following interpolations: 

l + 7y/4 + 3y 2 /2 + (3/2-7r/3)y 3 



BzppAv) 



1 + 3y/4 + (3/4 - n/6)y 2 



1 + lly/12 + (1/(6tt) + 1/6 + tt/18)?/ 2 + y 3 /(37r) + y 4 /(127r) 

= -(H) 

Notice that these choices are not unique. We found that the above expressions are 
ones of the simplest fitting functions that can reproduce the exact numerical results 
for the whole range of y, as we will see below. 

Figure [T] compares various approximations against the exact boost factor. In 
FigfJJa) (FigJU[b)), the exact boost factor B s{p) , Eq.© (Eq.flHJ) (solid line), its Tay- 
lor series up to the quadratic order 5 Tay i ors ( p ), Eq.® (dashed), 1/v approximation 
By V)S (p), Eq. ffTOj) (dotted), and our approximation -B a pp,s(p), Eq. ffTTj) (+), are shown 
as a function of y = a^pnx. Notice that our approximation -B a pp,s(p) completely falls 
together with its exact results. For the case of s-wave (p-wave) annihilation, the ac- 
curacy of the Taylor expansion Eq.([9]) decreases down to 99% for y = 0.45 (y = 0.64); 
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1 2 3 4 1 2 3 4 

y = a (n x) 1 ' 2 y = a (n x) 1 ' 2 

(a) (b) 

Figure 1: Exact boost factor B s f p \ , Eq.JT]) (Eq.©) (solid line), its Taylor series up to the quadratic 
order BTayior.s(p), Eq.(|9l) (dashed), 1/v approximation -Bi/u^p), Eq. lfTOl) (dotted), and our approxi- 
mation -B apPiS (p), Eq. ([TT|) (+), as a function of y = a\pKX. Our approximation B a pp, s(p) completely 
falls together with its exact results B s ( p y 



the range for the 1/v approximation Eq. ffTO]) to work at this level is y > 3.0 (y > 3.6). 
The in-between range 0.45 < y < 3.0 (0.64 < y < 3.6), where neither of the known 
approximations works, corresponds to 0.06 < a < 0.4 (0.08 < a < 0.5) for the typ- 
ical WIMP decoupling temperature T ~ m x /20. On the other hand, our ansatz 
Eq. (flT|) always reproduces the exact results with accuracy of less than 0.3% (0.9%). 



3 Relic Abundance with Sommerfeld Enhancement 

Let us discuss the computation of the relic abundance for the case where the anni- 
hilation cross section is enhanced by the Sommerfeld factor Eq.(T4]). As an analytic 
formalism for the computation of the relic abundance, we follow the standard freeze- 
out picture [U [5], with appropriate modifications [3TJ [32j [33J |3U [35]. Here, we show 
that the standard formalism works even for the case of Sommerfeld enhancement. 

The relic density of thermal relic particles x whose single production and decay 
are forbidden by some symmetry is determined by solving the Boltzmann equation, 

^ + 3^ x = -M(^-n^ eq ), (12) 

which describes the time evolution of the number density n x of the \ particles in the 
universe expanding at the rate H. Here n xm is the equilibrium value of n x , whose 
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non-relativistic limit is given by 

n x , eq = 9x (^) 3/2 ^ /T 5 (13) 

where g x denotes the internal degrees of freedom of the x particle. At high tem- 
peratures, the x particles are in thermal equilibrium. After T drops below m x , the 
number density n x exponentially decreases until the interaction rate T = n x (av) falls 
below the expansion rate in the radiation-dominated epoch H = (irT 2 / M P1 ) \J /90 , 
where Mpi = 2.4 x 10 18 GeV is the reduced Planck mass, and g* the number of the 
relativistic degrees of freedom. Then, the x particles are no longer kept in thermal 
equilibrium and the comoving number density becomes fixed. 

We can express the Boltzmann equation ([121 in terms of the dimensionless quan- 
tities Y x = n x /s and x = m x /T, where s = (2tt 2 /45)g*T 3 is the entropy density. 
Assuming that the universe expands adiabatically, the Boltzmann equation can be 
rewritten 



m x M P1 ^^(y 2 - Y 2 ) . (14) 



Introducing the variable A = Y x — Y X)Cq , we obtain 

dA _ _ «y„ _ 4, A(2y ^ + A) _ (15) 



dx dx V90 x 2 

The solution can be analytically derived in two extreme regimes. At temperatures 
above the freeze-out temperature T F , the x particles are in thermal equilibrium, so 
that the deviation of Y x from its equilibrium value K^eq is very small. Ignoring A 2 
and dA/dx, the solution is given by 

A~ _ X , (16) 

(8tt/ V90)m x M m y/gl(<jv) 

where we have used dY Xfiq /dx ~ — Y xeq for x 3> 1. Freeze-out occurs when Y x 



deviates from Y Xfiq : 



A(x F ) = cY x , eq (x F ) , (17) 

where xp = m x /T F , and c is a numerical constant of order of unity. Using the early 
time solution, Eq. (fl6|) . we obtain the recursive equation for determining the value of 

x F , 




,i/45 (av) 
x F = In ( \j —cm x Mpig x 



xg* 



(18) 

X=Xp 



For the standard s- and p-wave annihilation cross sections, the choice of c = \/z — 1 
is known to give a good agreement with the numerical results. We will see that this 
choice is still valid even with Sommerfeld enhancement. At temperatures below T F , 
the production term Y xm in the Boltzmann equation can be ignored. Therefore, the 
final relic abundance is found to be 

YJx oo) = =— , (19) 

X (4n/V90)m x M Pl I(x F ) 

where the annihilation integral is defined by 

/(»,) = fcte^M (20) 

For convenience, we express the final abundance in terms of Q x h 2 = m x soY x (x — > 
oo)h 2 1 p cr it, where s = 2.9 x 10 3 cm -3 is the present entropy density, and p crit = 
SM^H 2 is the critical density. The conversion factor of the annihilation integral to 
the relic density is given by 

n h 2 = 8 - 5 x 10-11 (2i) 

x I(x F ) GeV 2 ' 1 ' 

In sharp contrast to the exact boost factors, Eqs.([7j) and (jSJ), our ansatz B app , 
Eq. fllip . leads to simple analytic annihilation integrals. For s-wave annihilation, 

I(x F ) 




§a^f¥xp + 12 
(9 — 27i)Tra 2 XF , 

x 2(9-27r)a0fxF" + 9 
v/3(117 - 327r) 



(22) 



For p-wave annihilation, 



I(x F ) _ f°°j_B apPjP 



(>/\/7/7 J Xp dX x 



1 4av^F + (3 + vr 2 )a 2 + (21-7r 2 ) v ^Fa 



2x2 F 9 X 3 J 2 l8x F 36y/xF~ 

+ <3+ ; >W + 4 X 
144 V ay/Wx^ J 

In FigJ21 we show the normalized inverse freeze-out temperature function 
of the coupling constant a for the Sommerfeld-enhanced s-wave annihilation with 
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Figure 2: Normalized inverse freeze-out temperature as a function of a for the Sommerfeld- 
enhanced s-wave annihilation with a = 1.5 x 10~ 26 cm 3 /sec (a), and p-wave annihilation with 
b = 1.0 x 10~ 25 cm 3 /sec (b). Here we take m x = 100 GeV, g x = 2 and = 90. 



a = 1.5 x 10~ 26 cm 3 /sec (a), and p-wave annihilation with b = 1.0 x 10 -25 cm 3 /sec 
(b). Here we take m x = 100 GeV, g x = 2 and g* = 90. The plots illustrate that the 
freeze-out temperature decreases significantly as the coupling constant a increases. 

Figure [3]^a) compares the analytically estimated relic abundance using Eq. (122"]) 
(+) with the numerically computed result (solid line) for the same parameter set as 
in Figj2]^a). For the p-wave case, the approximated relic abundance using Eq. (l23|) 
(+) and its exact result (solid) are displayed in Figl^b). We also show the relic 
abundances for B s = 1 (B p = 1) (dashed) and those for the 1 jv -enhancement 
case (dotted). We emphasize that our analytic result successfully describes the 
transition from the non-enhanced case, B s — 1 (B p = 1) , to 1/v -enhancement, 
and overlaps with the numerically computed results even in the intermediate region 

3 x 10 -3 < a < 0.1, where neither of the limiting approximations works. In other 
words, the standard analytic method for the final abundance is also applicable to 
the Sommerfeld-enhanced annihilation cross section, which exhibits a non-trivial ve- 
locity dependence. The deviation of the analytically derived relic abundance from 
the numerical result is found to be always less than 1%. 

4 Conclusions 

In this paper, we have proposed an approximate analytic function for evaluating the 
boost factor caused by either s- or p-wave Sommerfeld-enhanced annihilation, Eq.flSJ), 
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a = 1.5 x 1(T 26 cm 3 /sec 




10 u 



p-wave 




(a) (b) 

Figure 3: Comparison of the analytically estimated relic abundance (+) with the numerically 
computed result (solid line). Here we take the same parameter set as in Figj2j The cases of B s = 1 
(B p = 1) (dashed) and 1/w-enhancement (dotted) are also shown. Notice that the analytically 
computed relic abundance using Eq. d^)) or (p?5|) (+) completely falls together with the numerical 
result (solid). 



assuming massless force carriers. For s-wave (p-wave) annihilation, our ansatz for 
the boost factor smoothly connects the limiting cases of B s = 1 (B p = 1) and 1/v- 
enhancement, and reproduces the exact values with accuracy of less than 0.3% (0.9%) 
in the whole range of the coupling constant a. We have applied the approximate 
boost factor -B apPiS (-B ap p,p) to the computation of the thermal relic abundance of 
particles whose s-wave (p-wave) annihilation is Sommerfeld-enhanced in the early 
universe, and derived a totally analytic formula for the final relic abundance in terms 
of the freeze-out temperature xf- Our findings show that the standard analytic 
method for the relic abundance attains accuracy of less than 1%. Our results are 
rather generic, and applicable not only to dark matter candidates but also to other 
particles that decouple non-relativistically from the thermal background in the early 
universe. 

Finally, we mention that chemical recoupling after kinetic decoupling could affect 
the relic abundance [21] • A more detailed discussion including the evolution after 
kinetic decoupling will appear elsewhere. 

After completion of this work, we received a preprint [37] that deals with a similar 
subject. 
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